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We investigate a bosonic quantum gas consisting of two interacting species in an optical lattice at 
zero and finite temperature. The equilibrium properties and dynamics of this system are obtained 
by means of the Gutzwiller mean-field method. In particular we model recent experiments where 
the ramp-up of the optical lattice occurs on a time scale comparable to the tunneling time of the 
bosons. We demonstrate the violation of adiabaticity of this process with respect to the many-body 
quantum states, and reproduce and explain the oscillations of the visibility as a function of ramp-up 
time, as seen in experiments. 
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I. INTRODUCTION 

Ultracold gases provide a powerful system for the ex- 
perimental investigation of interacting quantum many- 
body systems. In combination with optical lattices and 
tunable inter- and intraspecies interactions, degenerate 
quantum gases give insight into both strongly and weakly 
correlated regimes, and possible phase transitions. They 
allow the direct probing of fascinating phenomena like the 
superfluid to Mott-insulator (SF-MI) transition. In the 
breakthrough experiment [l[ this transition was shown 
for the first time, realizing predictions for the well-known 
Bose-Hubbard model @, Q in the case of a single bosonic 
species. 

Adding a second atomic species results in a wealth of 
quantum phases, clearly demonstrating the complexity 
of correlated ensembles. Currcntly,several experimental 
groups are working on Fermi-Bose , Fermi- Fermi 0- 
llOU and Bosc-Bosc mixtures in optical lattices as 

they are promising devices for studies of disorder flij ]. 
dipolar molecule formation [l5[ , and spin arrays [l6| . In 
Bose- Fermi mixtures [ljiEII and Bose-Bose mixtures [l9|, 
[20j . a supersolid phase has been predicted to exist. 

Recently the Florence group has realized a Bose-Bose 
mixture of 87 Rb and 41 K trapped in a 3D optical lattice 
[lH . They showed that the SF-MI transition of Rubidium 
is shifted towards shallower optical lattices when Potas- 
sium is present in the system. The same reduction of the 
visibility was seen in experiments where fermions were 
added to bosons in optical lattices 

These experiments are important steps towards the 
study of low-temperature properties of atomic mixtures. 
The experimental investigation of quantum many-body 
ground-states, for instance with spin ordering, requires 
the achievement of rather low temperatures in the lat- 
tice. For fermions this already poses problems before 
ramping up the optical lattice, since evaporative cooling 
becomes inefficient in the degenerate limit. Bosons, on 
the other hand, can be cooled to very low temperatures 
[2T| . However, ramping up the optical lattice can easily 
increase the temperature, or more precisely the entropy, 



again. This poses the question we address in this paper: 
under which conditions is the ramping up of the optical 
lattice sufficiently slow, such that the process is adia- 
batic. This question has been investigated before for the 
case of a single species bosonic gas; both by the mean- 
field technique we employ here and by methods tailored 
to the one-dimensional systems [23 - l24l | . 

Here we systematically map out the visibility of a Bose- 
Bose mixture as a function of ramp-up time. We repro- 
duce and explain the experimentally observed oscillations 
and relate them to the issue of adiabaticity. We take 
into account that this system breaks the rotational sym- 
metry of the trap, because the center of masses of the 
two species arc shifted with respect to each other due 
to the gravitational sag. We consider a two-dimensional 
system. In this way we keep the numerical effect man- 
ageable, and yet avoid the peculiarities associated with 
a one-dimensional system, like the absence of long range 
order even at zero temperature, which makes the extrap- 
olation of one-dimensional results to higher dimensions 
problematic. 

The timescalc for ramping the optical lattice in the 
experiment is usually determined such that the system 
ends up in the lowest band of the optical lattice, i.e. the 
ramping time is chosen large with respect to the band 
gap. Since the band gap is small for shallow optical lat- 
tices, sometimes an exponential ramp-up profile is cho- 
sen. However, this docs not guarantee that the ramping 
process is also adiabatic with respect to the many-body 
states in the lowest band [22|, [24[ • This is indeed not true 
and can be easily seen from the fact that the tunneling 
time for the atoms th op (s) = f j R 1( s ) * s 01 tne same order 
as the ramp-up time in most experiments. 

The results of our numerical simulations are in agree- 
ment with this qualitative argument. By simulating the 
ramp- up of the optical lattice using Gutzwiller mean-field 
theory, we observe that in the regime of deep optical lat- 
tices the ramp-up dynamics generally does not lead to 
the ground state of the system. The ground state of 
the system at those lattice depths typically contains a 
Mott insulating plateau with integer filling in the center. 
The states obtained after ramping up the lattice indeed 
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FIG. 1: (Color online) Modeling of the gravitational sag by 
two spatially separated parabolic trapping potentials. The 
lattice potential with the superimposed harmonic trapping 
potential (here drawn in one dimension only) is represented 
by a sinusoidal blue and red line. Different colors represent 
the two different bosonic species. 



lose the long-range superfluid order in this region. How- 
ever, unlike the Mi-phase, the particle densities are non- 
integer. We show that the exponential ramping profile is 
even less adiabatic than a linear increase of the depth of 
the optical lattice. 



By systematically tracing the visibility as a function 
of the ramp-up time, we find characteristic oscillations, 
which also have been observed in experiments with a 
single bosonic species [25| and a two-component bosonic 
mixture [ll|. The fact that we reproduce these oscilla- 
tions within the mean- field dynamics, where heating due 
to three-body collisions is not included, indicates that 
they are part of the real many-body dynamics. We in- 
terpret these oscillations as a competition between two 
effects. The first effect is the finite time the particles 
need to localize. For fast ramps, the number of dclocal- 
ized particles is far higher than the equilibrium value. 
However, phase coherence is generally not present. The 
second effect is a coupling to collective excitations, in 
particular the breathing mode of the system. Whereas 
the first effect enhances the visibility, the second effect 
destroys it, because the induced current leads to occupa- 
tion of finite momentum modes. The criterion for adi- 
abaticity, we deduce from this, is the saturation of the 
visibility as a function of ramp-up time. In this regime 
the superfluid fraction reaches its equilibrium value and 
the ramp-up is sufficiently slow such that no collective 
excitations are excited. The maximum of the visibility, 
as sometimes used as an experimental criterion, turns out 
to be not a good indication of adiabaticity. 



This paper is organized as follows: in section |H] wc 
present the two-component Bose-Hubbard model we in- 
vestigate and explain the Gutzwiller mean-field method 
we apply. In section IIIII we present results for T = 0. 
In section HVI we discuss effects of finite temperature. In 
section [V] we present conclusions and implications for 
the experiments. 



II. MODEL 

The two-species bosonic system we consider is de- 
scribed by the Bose-Hubbard Hamiltonian: 
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where the operator b\ a (b ia ) creates (annihilates) a bo- 
son of flavor a, a = a, b at site i. Since our goal 
is to model the experiment [Tl|, we choose the two 
species as 87 Rb and 41 K. The separation of the center 
of mass of the two clouds due to the gravitational sag 
is provided by two spatially separated parabolic traps 
V*™ P = \m a uj 1 a a^\xi — r Q | 2 , where d is the lattice con- 
stant, uj a the trap frequency, are the coordinates of the 
given lattice site i and r Q are the centers of the harmonic 
potentials, see Fig. [1] The frequencies of the parabolic 
trap are ljri, = 2tt x 36Hz and ojk = 2tt x 53Hz. The dis- 
tance between the two trap centers was chosen in such a 
way that the atomic clouds overlap only on a few lattice 
sites, as in the experiment ll|. The chemical potentials 



fx a are adjusted such that the resulting particle numbers 
correspond to the particle number ratio of the experiment 
N K /N Rb « 0.1 [ij (N K ps 30, N Rb ps 300). The parame- 
ters U a , U a bi J a indicate the intra-/interspecies Coulomb 
repulsion and the hopping amplitude. As input we use 
the experimental values of the Florence experiments ll| : 
lattice laser wavelength A^ = 1064 nm, scattering length 
dRb-K = 163ao, aiu, = 99ao, <xjf = 65ao, with ao the 
Bohr radius. The hopping constants J a and the inter- 
action parameters U a , U a b are calculated according to 
M: 
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with the dimensionlcss lattice depth s a = V"/E", where 
V" is the depth of the laser-induced potential and E" 
the recoil energy of species a. From the hopping con- 
stant we directly obtain the tunneling time of Rb-atoms 
ihop(s) = § jj^(s) ' w hich we plot as a function of lat- 
tice depth in Fig. [2] For the given wavelength the depth 
of the laser induced lattice potential is species-specific: 
Vl h = l.lVjf and Ef = 2.1E™. This results in 
SRb = 2.3sjf. In the following we use the short-hand 
notation s = s R b- Here wc neglect the influence of the 
parabolic confinement on the parameters. This is justi- 
fied because experimentally the potential energy differ- 
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FIG. 2: The tunneling time as a function of the lattice depth. 

ence between neighboring lattice wells is much smaller 
than the barrier height [27j . 

A. Gutzwiller Mean-Field approximation 

The Hamiltonian Eq. ([1]) can in principle be analyzed 
by Quantum Monte Carlo simulations [HI, [2i|. However, 
these simulations are restricted to equilibrium situations 
and cannot describe dynamical processes like the ramp- 
ing up of the optical lattice we consider here. Alterna- 
tively one can apply the time-evolving block decimation 
(TEBD) method to simulate the dynamics [24j . However, 
TEBD is restricted to one dimension. Here we simulate 
the experimental ramp-up dynamics within a mean-field 
approximation, which provides good qualitative results 
|22j, requiring a much lower computational effort than 
exact numerical alternatives. 

Within the Gutzwiller method a mean-field approxi- 
mation is applied to the operators b\ a (b ia ) in the hop- 
ping part of the Hamiltonian Eq. (p}. This leads to a 
decomposition of the lattice Hamiltonian into a sum over 
decoupled single-site Hamiltonians, which can be solved 
numerically: 

a—a.b 

-(fJ>a - V-™ p )h ia ] + U ab n ia n ib , (5) 

with ipi a = J2jn n ifija) i where the mean- field parame- 
ters (bj a ) (superfluid order parameters) have to be found 
selfconsistently. Since the Hamiltonian is a sum over 
on-site Hamiltonians, the many-body wave function is a 
product wave function over the lattice sites. However, the 
different sites are coupled by the superfluid order param- 
eter <f>i a . When the bosons are superfluid, the superfluid 
order parameter is nonzero and the phase is constant in 
space. This establishes superfluid long range order in the 
system. 

This mean-field approximation is exact in the weak- 
coupling limit [3^]. Moreover, it is exact in the limit of 
infinite dimensions. Corrections to the mean-field results 
scale like 1/z, where z is the number of neighbors [l9|, 
[3~i| . This means that on the cubic lattice, where z = 6, 



the mean-field theory is well-controlled. It provides still 
qualitatively good results in two dimensions (z = 4). 



B. Dynamics 

The Gutzwiller mean-field method can straightfor- 
wardly be generalized to time-dependent calculations 
0. At zero temperature, the total wave function 
is assumed to be a product wave function over the 
lattice sites: |*(t)) = ^ with = 

Each of the on-site wave-functions is evolved in time 
according to the local Schrodinger equation 

This constitutes a set of coupled non-linear differential 
equations for the /3i,n a ,n b (i), which have to be solved. 
The resulting mean-field dynamics conserves the average 
total particle numbers (N a ). In order to achieve number 
conservation numerically, the wave function on each lat- 
tice site is updated sequentially taking into account the 
already updated superfluid order parameters of the pre- 
vious sites. This method keeps the total particle number 
constant (deviations < 0.1 %o) 1321 , without need of a 
time-evolving chemical potential |22|. 



C. Finite temperature 

The Gutzwiller mean-field method is also readily gen- 
eralized to nonzero temperature. In this case the to- 
tal system is described by a density matrix p(t) . Ex- 
pectation values are obtained as (0)(t) = Tr[p(t)0]. 
In mean-field p(t) factorizes over the lattice sites as 

P(t) = UiPi(t), With pi(t) = E„ £ ^l^(t)>ii(V'nW| ) 

where \ip n (t))i are the local eigenstates of the Hamilto- 
nian 1hL 1 j AF {t) and Zi the on-site partition functions. This 
representation is also valid out of equilibrium. The local 
density matrix pi(t) obeys the Von Neumann equation 

ihd t p i (t) = [Ht IF (t),p i (t)} . (6) 

The von Neumann time evolution is unitary and hence 
preserves the weights in the density matrix. Thus the 

thermal weights in the density matrix — — in the ramp 
simulation are time-independent, where E n are the initial 
eigenvalues of the system and Zi is the initial partition 
function. The initial inverse temperature (3 = 1/HbT is 
set to an experimentally reasonable value. Since in the 
Gutzwiller approach the density matrix factorizes over 
the lattice sites, this puts severe restrictions on the pos- 
sibility of (local) thcrmalization. In order to account for 
dissipation, such as induced by three-body losses, one has 
to apply the Lindblad equation [33[. 
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Mean-field theory predicts the existence of long range 
order at low but nonzero temperatures in any spatial di- 
mension. This means that the correct physics in two- 
dimensions is not fully recovered at finite temperatures. 
In three dimensions, mean-field theory describes the cor- 
rect behavior of the system at zero and nonzero tem- 
perature. Our two-dimensional simulations can thus be 
viewed as a qualitative description of the experiments in 

M. 



III. RESULTS FOR T — 

At low lattice depth the system does not fulfill the 
single-band and tight-binding approximation required for 
the Bose-Hubbard Hamiltonian in Eq. ([1]) to be valid. 
We therefore start our calculations at s = 5 where 
87 Rb is already far in the tight-binding and the lowest 
band regime. The corresponding lattice depth for 41 K is 
sk = 2.17. Since the ratio between the bandwidth 4Jr- 



and the bandgap E g 



2J~SkE r is small for the initial 



sk (4:Jk I Egap ~ 0.03), the single-band approximation 
is also satisfied for 41 K. We assume that the wave func- 
tion still corresponds to ground state when the ramp in 
the experiment reaches this initial lattice depth. This is 
a valid assumption, since interactions are still weak for 
s = 5 and the ramping time in the experiments is chosen 
adiabatic with respect to the band gap, which guarantees 
that the particles remain in the lowest band. Starting 
with s = 4 and s = 6 indeed did not change our results. 

Although the ramp of the lattice affects both species, 
only effects of the ramp-up dynamics on the 87 Rb-atoms 
are discussed in the following, like in the experiments 

Da- 

To investigate the adiabaticity of the lattice-ramp, we 
now first present results for the density distribution in 
real-space as well as in momentum space. Later on we 
focus on the visibility. 



A. Density profiles 

1. Real space 

To investigate the effect of the ramping dynamics we 
first compare the density distributions in real-space af- 
ter the lattice ramp-up with static density profiles at the 
corresponding final lattice depth. The static profiles orig- 
inate from the ground state at T = for a given s. 

For small final s (s < 10) the density profiles agree 
perfectly with the static ones independently of the ramp- 
ing time. For larger s, longer ramping times are needed 
to achieve good agreement between the profiles. In par- 
ticular, for s > 15 in equilibrium a Mott plateau ap- 
pears in the center of the trap. However the dynami- 
cally evolved wave function remains a superposition of 
various Fock states for fast ramps. This can be under- 
stood by the following argument. In the limit of an in- 
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FIG. 3: (Color online) Particle density and superfiuid order 
parameter of 87 Rb along a cut in ^-direction through the cen- 
ter of the trap (y — 0). (a) For s = 15 and ramp duration 
t = 50 ms (red dashed line) the particles in the center of the 
trap are in the "frozen phase" in contrast to the ground state 
profile (blue solid line), where a Mott-plateau is present. This 
is indicated by the noninteger density n Rb and nonzero local 
superfiuid order parameter |(6r&)| (c). (b) For s = 22 and 
t — 50 (red dashed line) a density wave is induced due to 
the fast ramp and high nonadiabaticity. Global parameters: 
L = 45, N Rb = 303, N K = 30, Um,-k = 1.93*7*6. Parame- 
ters for (a) and (c): J Rb = 0. 02URt, Parameters for (b) and 
(d): J Rb = 0.00417™. 



stantaneous increase of the lattice depth the superfiuid 
wave function consisting of local superpositions of var- 
ious Fock states remains unchanged. Since the hop- 
ping J decreases exponentially with increasing lattice 
depth, after this sudden step to a deep lattice the time- 
evolution operator consists mainly out of the interaction 
part Uh{n — l)/2 which rotates the phase of each Fock 
state independently, but leaves the local wavefunction in 
a superposition of multiple Fock states. The situation 
is similar for fast ramp and leads to noninteger particle 
density (see Fig. [3ja)(b)) and nonzero local superfiuid 
order parameter (see Fig. [3jc)(d)). Although the abso- 
lute value of the superfiuid parameter is not vanishing, 
long-range order is destroyed as the complex phases of 
(pi) are not constant over the lattice any more. Cor- 
relations between the phases, however, can be partially 
recovered for certain ramp-up times. This corresponds 



to the collapse-and- revival physics [251 ]: after time inter- 
vals of length 2irh/U all the individual phases are back 
in phase and global phase coherence is restored. Away 
from the revival times, due to the vanishing coherence 
this phase is not a SF-phase, we refer to it as a "frozen" 
phase. 

Ramping the lattice to higher s with the same ramp-up 
time means effectively faster increase of the lattice depth 
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FIG. 4: (Color online) Particle density in momentum space 
along a cut in (l,0)-direction for k y — 0. (a) In the static 
calculations (blue solid line) the central peak at k = (0, 0) is 
present, due to the SF shell, and small side peaks at nonzero k 
due to the Mi-core in the middle of the trap. Fast ramps with 
t = 50 ms (red dashed line) drive the particles into a frozen 
phase with lowered phase coherence, which is mirrored in a 
decreased central peak, (b) Ramping within the same time 
to a deeper lattice is highly nonadiabatic. The central peak 
becomes smaller than the side peaks or even the background 
value n(k = (n, 7r)) and the ±k symmetry disappears. 



per time unit. This causes not only significant deviations 
from the static density profiles but also the formation of 
density waves (Fig. EJb)). This behavior is observed for 
ramp-up times below 100 ms. The density waves are a 
further evidence that the final state reached after the 
lattice ramp is not necessarily the ground state. 

We finally examine the rotational symmetry of the den- 
sity profile of the Rb atoms. Due to the presence of the K 
species (and hence additional repulsive interactions), the 
rotational symmetry of the Rb cloud is always broken. 
However, this effect is only very pronounced for large 
s. Deep in the MI regime the compressible supcrnuid 
boundary layer, adjacent to the cloud of K-atoms, disap- 
pears because of the strong repulsive force experienced 
by the K atoms. The corresponding superfluid order pa- 
rameter vanishes in the ground state. This is different in 
the case of the time-evolved profiles: as the 87 Rb atoms 
for t = 50 s = 22 ms are in the "frozen" phase, the break- 
ing of circular symmetry is far less pronounced as in the 
static case, where the system is already Mott-insulating. 



2. Momentum space 

More information regarding superfluid long-range or- 
der is available in momentum space. Besides, in contrast 
to the real-space particle distribution, information on the 
momentum distribution is experimentally well accessible 
by time-of-flight measurements. The nonadiabaticity and 
oscillations in the density profile are mirrored in the mo- 
mentum distribution. To accentuate the global behavior 
as a function of the ramp-up time and the final lattice 
depth we show only the values averaged over an equili- 



bration time of 16 ms. 

Whereas for a homogeneous system of nonintcracting 
particles the momentum distribution is a single delta- 
peak at k = (0,0), the momentum distribution of a 
trapped system of interacting particles has a constant 
background due to non-condensed particles and shows 
broadening because of the trap. This latter effect dom- 
inates over the broadening by short-range fluctuations, 
which are neglected in the Gutzwiller approximation. 

As pointed out in the previous subsection, for short 
ramp-up times and deep lattices the local superfluid order 
parameter remains finite in the region where the static 
calculation predicts a Mi-plateau. Particles stay "frozen" 
in a superfluid-like phase indicated by noninteger parti- 
cle densities and finite local superfluid order parameter. 
Although in the dynamic case all Rb atoms seem to be 
superfluid, the phase coherence is lost. This is shown in 
Fig.[H The reduced central peak in the n(fc)-profile com- 
pared to the static data clearly indicates destroyed long- 
range order in the system (Fig. Ufa)). In extremely nona- 
diabatic cases the central peak becomes smaller than the 
side peaks or even the background value n(k = (tt,tt)). 
This leads to a significant broadening of the momentum 
distribution. Additionally the ±fc symmetry breaks down 
(Fig. lib)). 



B. Visibility 

The most convenient way to compare the experimental 
momentum distribution with our theoretical results is to 
calculate the visibility rj for the 87 Rb atoms as a function 
of lattice depth s 



n s (fc = (0,0))-n s (fc = (7T,7r)) 
n s (k = (0,0)) + n s (k = (rr,7r)) 



(7) 



Here n s (k) corresponds to the spatial Fourier transform 
of (b\bj)(s). In the experimental procedure [ll[ the 
height of the first order peaks is compared with the min- 
imum in n(k) at the same distance from the central zero 
order peak, to divide out the contribution of the Wan- 
nier function. Our calculations are performed within the 
tight-binding model Eq. (J]) and do not include the shape 
of the Wannier function. We therefore calculate the vis- 
ibility by comparing the central peak with the minimum 
at the edge of the Brillouin-zone. 

The n(k) values are extracted in the experiment (TlT | by 
integrating finite square areas around the peaks instead 
of taking single values. Application of this method to 
our theoretical data shifts the visibility for all ramp-up 
times. This is a small quantitative effect and depends on 
the extent of the integration area. It has no effect on the 
conclusions drawn in this paper. 

In order to investigate the reduction of the visibility in 
a systematic way, we now subsequently analyze the role 
of the second species, the ramp-up profile of the lattice 
and the ramp-up time. 
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FIG. 5: (Color online) Visibility for pure 87 Rb (blue solid line) 
and in the presence of 41 K (red dots) after a. t — 300 ms ramp 
(a) and in the static case (b). In both cases the visibility of 
pure 87 Rb is higher. 



1. Effect of the second species 

It is observed experimentally that addition of 41 K to a 
system of 87 Rb-particles reduces the phase coherence and 
the visibility [ll| . Our simulations reproduce this behav- 
ior for similar parameters and particle number ratios as 
in the experiment, i.e. Nrj, ps 303 and Nk ~ 30. We ob- 
served this effect for all simulated ramping times (from 50 
to 300 ms with a linear ramp-up profile) as well as in the 
static case (Fig. [5] (a), (b)). In particular, the destructive 
effect of the 41 K on the phase correlations of the super- 
fluid order parameter of 87 Rb is more pronounced in the 
dynamic case than in the static one. This indicates that 
the second species enhances the noncquilibrium induced 
by the lattice ramp. 

The reason for the lower 87 Rb visibility in the presence 
of 41 K in the static case is the following. The repulsive 
interaction between the species pushes the 87 Rb atoms 
out of the overlap region and enhances the 87 Rb density 
in the trap center. This increases the interaction energy 
of the 87 Rb and brings the particles closer to the SF-MI 
transition, thus decreasing the coherence of the system 
in the static case [23| . 

The dynamic lattice ramp to deep lattices leads even 
to a lower visibility than in the static case. This can be 
explained by the following argument. As we will show 
below, one of the reasons for the nonadiabaticity of the 
ramp-up is the excitation of collective modes. This effect 
results from the interplay between the intraspecies repul- 
sion and the increase of the lattice depth, which squeezes 
the atomic cloud. When 41 K is present, the clouds exert 
a repulsive force on each other, even though the over- 
lap region is small. This leads to additional collective 
modes in the system, which lowers the phase coherence 
and hence the visibility. Below we will quantify this state- 
ment. This mechanism could explain why in experiments 
always a reduced visibility is seen when a second species 
is added. 



FIG. 6: (Color online) Effect of the ramping profile on the 
visibility of 87 Rb in the presence of 41 K. (a) The resulting 
visibility after t — 75 ms ramp with the linear ramp (blue 
solid line) lies closer to the adiabatic static visibility (red solid 
line) than the visibility after ramp with the exponential profile 
(green dashed line), (b) Results with the same profiles but 
for t = 200 ms ramp. 

2. Effect of the ramp-up profile 

In this subsection we investigate the effect of the time 
profile of the lattice ramp on adiabaticity. Our motiva- 
tion is that in experiments the ramping profiles are usu- 
ally of exponential shape to keep particles in the lowest 
band at the beginning of the ramp. For a better com- 
parison between the simulated linear profile and the one 
used in experiment we also used the same shape as in 

s'(t') = (e^-l)^— , (8) 
eo.4 — l 

where s is the final lattice depth and t the ramping 
time. We performed calculations for 50to 300 ms ramp- 
up times. These simulations demonstrate that the expo- 
nential ramp leads to a lower visibility than the linear 
profile (Fig. |6]). The highly nonadiabatic ramp- up time 
of 75 ms with exponential profile leads even to negative 
visibilities (Fig.UJa)). This indicates that the system dis- 
plays collapse-and-revival physics for these short ramping 
times. Extending the ramping time to t = 200 ms results 
in an almost adiabatic linear ramp for s < 15 as the 
calculated visibility corresponds very well to the static 
(Fig. (Ha)). The visibility of 87 Rb after 200 ms exponen- 
tial ramp is higher than after 75 ms but still lower than 
the linear results. For deep lattices with s = 22 the expo- 
nential ramp-up profile leads to a higher visibility than 
the linear ramp. However, this is due to the fact that the 
exponential ramp-up is still highly non-adiabatic for this 
ramping time, leading to a large SF-fraction, whereas the 
linear ramp is closer to being adiabatic. 

The observation that the exponential ramp-up pro- 
file leads to less adiabatic behavior is explained by the 
fact that this profile has a low ramp-up velocity for the 
smaller lattice-depths, where the interaction only plays a 
minor role, whereas for the higher lattice depths, where 
many-body effects become important, the ramping ve- 
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FIG. 7: (Color online) Calculated visibilities for 125 ms (green 
squares) 100 ms (pink crosses) and 50 ms (blue dots) ramping 
durations compared to the static visibility (black dots, not 
connected by line) . The inset demonstrates that for s < 17 the 
visibility after a 125 ms ramp lies closer to the static visibility 
than after a 100 ms and vice versa for deeper lattices. The 
increase of the visibility for deep lattices for t — 50 ms is due 
to the highly non-adiabatic revival of the coherence. 



locity is very high. For this reason further investigations 
were performed exclusively with a linear ramp-up shape. 




300 



FIG. 8: (Color online) Visibility as a function of ramp-up 
time. For s = 7 (red solid line) all considered ramping times 
are adiabatic and the visibility is independent of the ramp- 
up time. For s = 14 (blue asterisks) longer ramp leads to 
the ground state visibility, indicating 300 ms as the appro- 
priate ramp-up time. For s > 17 oscillations are induced, as 
explained in the text below. 



3. Effect of the ramp-up time 

Finally we investigate the effect of the ramp-up time on 
the visibility. Whereas one would expect that a slower 
ramp automatically enhances the adiabaticity and the 
visibility, in the simulations we observe that only the dy- 
namic real-space particle distributions become similar to 
the static ones. The visibility, and therefore the phase 
coherence, remains different from the ground state one. 
Especially for intermediate ramping times the behavior is 
counterintuitive (Fig. [7]). While for lattice depth s < 15 
the visibility increases monotonically with the ramp-up 
time (Fig.[7]inset), deeper in the Mi-regime this tendency 
is washed out. This is surprising as the longer ramp is 
expected to be more adiabatic. Plotting the visibility as 
a function of ramp-up time for a fixed s demonstrates 
that the effect of the ramping time is not the same for all 
s (Fig. [5]). We distinguish two different regimes. Ramp- 
ing up to shallow lattices s < 12 (h/J(s = 12) w 5.4 ms) 
is adiabatic on all time scales and thus not affected by 
the ramp duration. In the regime of final lattice depth 
around the SF-MI transition and in Mi-phase, oscilla- 
tions in the visibility occur. 

These oscillations were also observed experimentally 
[ill . |25| but not yet explained. In the next subsection we 
give an explanation in terms of coupling to the collective 
modes of the system. Before this we add some remarks 
on the figures. The oscillating behavior depends on the 
ramp-up profile (Fig. IHa)(b)). The exponential shape 
shifts the oscillations to longer ramp-up times. The pres- 
ence of 41 K has only a minor effect on the oscillations: the 
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FIG. 9: (Color online) Effect of the ramping profile on the 
visibility oscillations (a), (b) and effect of the second species 
(c), (d). For intermediate lattices e.g. s = 16 the exponen- 
tial ramp-up profile (blue dashed line) damps the oscillations, 
which are present for the linear profile (red solid line) (a) . For 
s = 22 and exponential ramp (b) a non-adiabatic maximum 
appears. The presence of a second species (blue dashed line) 
(c) , (d) only induces a global shift of the visibility. 



second species only leads to a global shift in the visibil- 
ity without changing the position of maxima and minima 
(Fig. El[c)(d)). This is consistent with the experimental 
findings 
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FIG. 10: (Color online) (a) Oscillations of the radius of the 
Rb-cloud during the waiting time of 16 ms after a to = 125 ms 
ramp to s = 22. The data is fitted by f(t) = 0.135 sin(0.18t + 
5.19)exp-°' 031(t - to) +48.857 (blue line). The sinusoidal form 
corresponds to the excitation of the breathing mode, (b) The 
current at s — 22, averaged over the waiting time and renor- 
malized with respect to SF-fraction, induced after different 
ramp-up times. The maximum of the renormalized current 
indicates the regime with the maximal coupling to the col- 
lective modes. The second species leads to an enhanced cou- 
pling for all ramp-up times, (c) The SF-fraction averaged over 
the waiting time as a function of ramping time at s = 22. 
The "freezing" in the SF-phase is dominant for fast ramps. 
For longer ramp-up the system approaches the ground state 
(static) value and the ramp is getting more adiabatic. The 
41 K slightly reduces the SF-fraction. 



4- Explanation of visibility oscillations 

We now turn to the explanation of the oscillations in 
the visibility. A closer look at the simulations shows that 
the only component contributing to the visibility which 
oscillates in time at a fixed s is n(k = (0, 0)). The particle 
density in momentum space is given by 
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where L 2 is the total number of lattice sites, J"$(rj) and 
4>i(4>j) are the absolute value and the phase of the lo- 
cal superfluid order parameter on a site i (j) respec- 
tively: (b\) = r i e~ l ^' i . We can separate the contribu- 
tion of the absolute value and the phase to the visibil- 
ity. The absolute value r, shows monotonic behavior 
and is continuously decreasing if the ramp-up time is re- 
duced. This is measured by the Gutzwiller SF-fraction 
fc = T,i\{t>i)\ 2 / N - In Fi g- [Hid) the SF-fraction is con- 
tinuously decreasing and always higher than the equi- 
librium value. Only for long ramping times (t = 300 
ms), the dynamical and static SF-fraction approach each 
other. This means that for the short ramping times 
the visibility is dominated by this anomalously large SF- 
fraction. However, for very fast ramping times the SF- 
fraction is high, but the phases are uncorrelated, leading 



to a low visibility. With increasing ramping time the 
phase coherence builds up and compensates the decay of 
the SF-fraction leading to an increasingly higher visibil- 
ity. 

For even longer ramping times the visibility decreases 
again. We explain this by an enhanced coupling to the 
collective breathing mode of the system induced by the 
lattice ramp. Increasing the lattice depth results in an 
increased ratio of the strength of the harmonic trapping 
potential and the hopping constant J. This forces the 
particles to move towards the center of the parabolic trap 
and leads to a higher occupancy in the middle of the trap. 
At the same time the particles experience a higher repul- 
sion as U(s) grows with the lattice depth (see Eq. ((3J)). 
This repulsion acts against the increasing population and 
induces a reverse flow. The interplay between these two 
mechanisms yields the collective modes. 

We indeed find numerical evidence for these collective 
oscillations by observing the cloud size R 2 = (r 2 ) — (r) 2 
during the waiting time after the ramp. The sinusoidal 
oscillations of R 2 indicate the collective movement of par- 
ticles within the breathing mode in Fig. HPT a). The data 
for other ramp-up times can also be fitted similarly. This 
leads to the conclusion that only the breathing mode is 
excited. Depending on the ramping time the amplitude 
of the oscillation changes. In particular, the amplitude is 
continuously decreasing when the ramping time is made 
longer. However, this is mainly due to the decreasing SF- 
fraction, which reduces the number of mobile particles. 
Renormalizing the amplitude by the SF-fraction leads to 
a peak at the position of the minimum of the visibil- 
ity, which evidences that the coupling to the modes is 
responsible for this minimum. 

Analysis of the total current in the system I = 
Yltij) — bji>i)\ 2 further clarifies this. The total cur- 

rent is decaying because of the decaying SF-fraction. 
To investigate the relative motion of the mobile parti- 
cles we therefore renormalize the total current by f 2 : 
I T = I /N f 2 . This function again shows a clear maximum 
at the position of the visibility minimum (see Fig. llOf c)). 
As the coupling to collective modes destroys the phase co- 
herence, this fully agrees with the minimum in visibility 
for a 125 ms ramp, see Fig. Eld). 

For t > 150 ms the SF-fraction approximates the static 
value and the collective modes are less excited. This is 
the most adiabatic ramping regime. 

For the exponential ramp-up profile, the oscillations 
are shifted (Fig.[S£a)(b)). This is because the SF-fraction 
remains anomalously high even for long ramping times 
and the visibility is dominated by this effect. In par- 
ticular the maximum at t ramp = 200 ms for s = 22 is 
explained by the high SF-fraction and is thus a highly 
non-adiabatic point. The decrease of the visibility at 
t = 300 ms for s = 22 is explained by an enhanced cou- 
pling to collective excitations. 

The fact that the presence of a second species leads 
only to a small shift in the visibility is explained by the 
observation that the overlap of the atomic clouds is very 



small and the modes are mainly excited by the increased 
repulsion between the 87 Rb particles when the optical 
lattice is ramped up. However, the second species induces 
additional modes, which lower the visibility. This is seen 
in the higher renormalized current in Fig. HOf b) for the 
mixture compared to the single-species system. 

The additional induced modes in the system due to the 
presence of K not only explain the lower visibility, but 
also offer an explanation for the experimental observa- 
tion that adding a second species leads to a broadening 
of the momentum profile beyond a certain lattice depth 
[TTI | . The presence of collective modes leads to macro- 
scopic occupation of single particle states with nonzero 
momentum and hence broadens the momentum distri- 
bution. It is worth noting that for this explanation the 
amount of spatial overlap of the two species is less im- 
portant: as long as the two clouds touch, they can exert 
a force on each other. This explains why this effect was 
already found for widely separated components. 

IV. RESULTS FOR NONZERO T 

In order to understand the effect of finite temper- 
ature in the experiments we also perform simulations 
for this case. We perform simulations for initial tem- 
peratures of 19 nK (fcsT 2 = 2.2Jub,s=5) and 12.6 nK 
(ksTi = 1.5Jb6,«=s). This is in the range of typical ex- 
perimental temperatures, which can be estimated as de- 
tailed in appendix [A] To investigate the adiabaticity of 
the ramp at finite temperature we again compare time- 
dependent ramp-up simulations with static results. The 
latter correspond to the ensemble in thermal equilibrium 
at a final lattice depth with effective inverse temperature 
$ = 1/ksT. The effective temperature is chosen such 
that the entropy of the static system equals the initial 
entropy of the ramped system. The static results thus 
represent an adiabatically ramped and completely ther- 
malized ensemble. Accordingly, the static density ma- 
trix is p = J2 n e Zi" \ E n)a{E n \, where E n , \E n ) t are the 
eigenvalues and eigenstates of the Hamiltonian at site i 
respectively. 

A. Density profiles 

We first investigate the density profiles in real-space 
at UbTi = l-5Jfli, :S= 5. Ramping in 50 ms to s — 12 pro- 
vides a density profile similar to the static thermalized re- 
sult with effective temperature fc^Ti = 0.95JRb. s= 5 (see 
Fig. [TTfa)). The local supcrfluid order parameter, how- 
ever, differs from the value in thermal equilibrium (see 
Fig. ITTT c)). In the center of the 87 Rb cloud the local su- 
perfluid order parameter is reduced. Hence the system 
after the ramp does not correspond to an adiabatically 
ramped and thermalized ensemble. At deeper lattices the 
dynamic density profiles do not fit the thermal static dis- 
tribution. Fig. Illf b) demonstrates the squeezing of the 
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FIG. 11: (Color online) Finite temperature results for the par- 
ticle density and superfluid order parameter of 87 Rb at final 
lattice depth s = 12, 16 along a cut in x-direction through 
the center of the trap (y = 0). All dynamic profiles are 
calculated for the temperature fcsTi = 1.5J R b,s=5- (a) For 
50 ms ramp (red dashed line) the results correspond well to 
the static thermalized ensemble (blue solid line) with effective 
temperature fcsTi = 0.95J Rb s=5 . However, the local super- 
fluid order parameter is reduced in the center of the atomic 
cloud (c). (b) After the ramp to s = 16 the dynamic profiles 
differ from the static ones with corresponding effective tem- 
perature fcsT2 = 0.42Jr6 iS =5. A density plateau is formed 
in the center of the trap at a noninteger density. For 50 
ms ramp (red solid line) density waves appear around this 
plateau, (d) The superfluid order parameter (red solid line) 
is peaked in the region where the density waves appear. Pa- 
rameters: L = 60, N Rb — 303, Nk = 30, Uwu-k = l-QSUm, 
Parameters for (a) and (c): Jut = 0. 24U R b, Parameters for 
(b) and (d): J Rb = 0.02U Rb . 

dynamic profile for s = 16 compared to the static one 
with effective temperature ^bTi — 0.42 Jnb.s=5- Similar 
as in the T = case, the fast 50 ms ramp is nonadia- 
batic and induces density waves around the plateau (see 
Fig. fTTT b)). The peaked local superfluid order parameter 
in Fig. lllf d) corresponds to the region where the density 
waves appear. For t = 300 ms in Fig. [TlT b) a plateau is 
formed at non-integer density in the center of the trap. 
At the same time the corresponding local superfluid or- 
der parameters vanish which indicates a formation of a 
normal phase instead of a Mi-plateau as found previously 
in the T = case. 



B. Visibility and Oscillations 

As shown in the previous subsection, with increasing 
temperature and lattice depth the condensate depletes. 
This also lowers the visibility. In Fig. Q2] this behavior 
is exemplified for a 300 ms ramp. For low temperature 
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in the center of the trap for slow lattice ramp compared 
to fast ramp (Fig. QTJd)), the phase coherence in the lat- 
ter case is almost completely destroyed. This leads to 
a higher visibility for a 300 ms ramp than for a 50 ms. 
For ksT2 = 2.2 Jnb,s=5 the visibility oscillations are sup- 
pressed and the minimum at t — 200 ms disappears. We 
can understand the disappearance of the minimum in the 
visibility in a qualitative way by comparing the excita- 
tion energy of the breathing mode with the temperature. 
We indeed find that the temperature here is higher than 
the excitation energy, meaning that the mode is already 
thermally occupied and that the ramp of the lattice has 
loss effect. 



FIG. 12: (Color online) Visibility after 300 ms ramp-up for 
k B Ti = 1.5 Jia,s=5 (green dots) and k B T 2 = 2.2J fl6jS=5 (blue 
triangles) compared with T — (red solid line) . The visibility 
decreases with increasing temperature. 
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FIG. 13: (Color online) Visibility as a function of ramp- 
up time for different temperatures and final lattice depths, 
(a) At shallow lattices for T = (red solid line) oscilla- 
tions appear with a visibility minimum at t — 200 ms. For 
fcsTi = 1.5J_H6 iS= 5 (blue dashed line) the global visibility 
and the oscillation amplitude is reduced. Higher tempera- 
ture (fcsT2 = 2.2Jr6 >s= 5, green dotted line) blurs the os- 
cillations, (b) For s — 19 the shape of the oscillations at 
T = is mainly conserved for k B T\ = 1.5J Rb a=5 . The graph 
is shifted to lower visibility and longer ramp-up times. At 
k B T 2 = 2.2jRb, a =5 only one maximum at t = 200 ms is 
present. 



ksTi = 1.5 Jm, s =5 and shallow lattices s < 12 the visibil- 
ity is hardly changed compared to T = 0. With increas- 
ing lattice depth, however, the superfluid order param- 
eter vanishes in the trap center, leading to a decreased 
visibility. For higher temperature (ksT 2 = 2.2Jm, tS=5 ) 
the fraction of the atoms in the normal phase increases, 
thus lowering the visibility further. 

The visibility oscillations are affected as well. When 
temperature and final lattice depth are sufficiently small, 
the only effect of the temperature is to reduce the visi- 
bility (Fig.lTJa) k B T x = 1.5Jk6, s=5 ). Although the ab- 
solute value of the local superfluid parameter is reduced 



V. DISCUSSION AND CONCLUSION 

Using a time-dependent Gutzwiller model for an inter- 
acting Bose-Bose mixture we investigated the ramp-up 
of the optical lattice for zero and at finite temperatures. 

The non-adiabaticity of the lattice ramp was ana- 
lyzed by comparing density profiles in real and momen- 
tum space and by studying the visibility. The adiabatic 
regime is reached when the density profiles as well as 
visibility agree with the equilibrium results. 

We have shown that a ramp-up of the optical lattice 
carried out on a time scale comparable to the tunneling 
time does not necessarily provide the ground state of the 
system. Depending on the ramp-up time, ramping the 
lattice at T — leads to the trapping of the particles in 
a "frozen" phase with non-integer particle number and 
nonzero local superfluid order parameter but vanishing 
global phase coherence. Ramping the lattice at finite 
temperature additionally causes a depletion of the con- 
densate. The latter grows with increasing temperature. 
Both ramp effects lead to decreased visibility. 

The ground state visibility was only reached for shal- 
low lattices within the investigated ramping times. The 
fact that one needs rather long times to be completely 
adiabatic for deep lattices, is rooted in the critical slow- 
ing down of the hopping at the SF-MI-phase boundary 
at T = or the SF-normal-phasc at finite temperature. 

For T = we demonstrated in addition that the linear 
lattice ramp is more adiabatic than the exponential. We 
found that a longer ramp-up time does not naturally lead 
to a better visibility. In fact, depending on the final 
lattice depth oscillations may occur for T = and low 
temperatures. 

We explain these oscillations by a coupling of the ramp- 
up process to the collective modes of the system. This 
is consistent with the appearance of density waves in 
the system. Lowered superfluidity and a larger normal 
phase prevents collective excitations at higher temper- 
atures. One of our main results is that the maximum 
in the visibility is not a good indication of adiabatic- 
ity. This regime is in fact highly non-adiabatic, since the 
maximum is caused by an anomalously large SF-fraction 
induced by the short ramping time. The region where 
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the visibility saturates is the most adiabatic. However 
in our approach three body collisions and heating is not 
included. These processes become relevant at long time 
scales and also lead to non-adiabaticity. 

In experiments, the presence of a second species de- 
stroys the phase coherence of the majority species, lead- 
ing to a decreased visibility. In contrary, previous theo- 
retical static calculations predict either enhanced or de- 
creased longrange order, depending on the actual par- 
ticle ratio J23(. Here we observed that for T = the 
dynamical ramp induces additional non-equilibrium and 
leads to a more pronounced visibility decay in the re- 
gion where static calculation also predicts lowering. We 
explained this in terms of an enhanced current in the 
system. This supplementary visibility reduction could 
be the reason for the experimental findings. 
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Appendix A: Estimation of the Temperature in a 
weak lattice 

The ramp-up in our work starts at s = 5 and not at 
s = 0, as in the experiment, due to the required tight- 
binding regime of the Hubbard Model. Therefore, start- 
ing from s = 5 we have to recalculate the temperature 
based on the initial experimental temperature 73 nK [ll| 
before the ramp. It can be shown by the following ar- 
gument that the ramp-up of the optical lattice cools the 
system down. For a first estimate let us assume that the 
lattice is ramped up adiabatically from s = to s = 5 
and that the particles are noninteracting. The initial 



slow ramp of the exponential or any other ramping pro- 
file realizes the first condition. The second condition is 
strictly satisfied only in shallow lattices. For s = 5, where 
Um ~ Jm, it can be assumed to be roughly satisfied and 
leads to the right temperature range. 

In the case of an adiabatic lattice ramp the entropy of 
the system S = -k B J2k y m [l + + n k) - nk In rife J 

remains constant. As the experiments start from s = 0, 
the initial dispersion corresponds to the free particle dis- 
persion e\ = h 2 k 2 /2m. The final situation is a tight- 
binding case and the particle density is evaluated with 
the dispersion e{ = — 2J(cos(k x a) + cos(k y a)), where 
a is the lattice constant. As the noninteracting con- 
densate is located around k = (0,0) the tight-binding 
dispersion can be approximated by a Taylor series as 
e{, w J(fca) 2 . The initial and final particle densities 
7i£ = (e /3£fc — 1) _1 have the same functional depen- 
dence on momentum k. The sum over functions with the 
same functional ^-dependence in entropies can only re- 
main constant when these functions are identical. This is 
realized when the particle density for the initial lattice 
depth corresponds to the final one. From the equality of 
initial and final particle densities the criterion for the fi- 
nal temperature follows /T l = st/s].. From here the 
final temperature can be estimated T f = T i J\ 2 m/2ti 2 , 
where A is the laser wave length and J is the hopping 
amplitude at the final lattice depth [34j . 

For a lattice ramp to s = 5 the initial temperature 
of 73 nK is lowered by a factor 2. As for this lattice 
depth the system is in the SF regime, the thermal energy 
competes with the hopping energy J. At this point the 
dimensionality of the experiment has to be taken into 
account. As the hopping scales with the number z of the 
next neighbors the ratio / zJ for our 2D system studied 
here should correspond to the ratio of the experimental 
3D system. This lowers the calculated temperature by an 
additional factor 2/3, providing the effective temperature 
T* = 24 nK for our equivalent 2D system. 
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